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Abstract. The quantum cosmological wavefunction for a quadratic gravity theory derived 
from the heterotic string effective action is obtained near the inflationary epoch and during 
the initial Planck era. Neglecting derivatives with respect to the scalar field, the wavefunction 
would satisfy a third-order differential equation near the inflationary epoch which has a solution 
that is singular in the scale factor limit a(t) — > 0. When scalar field derivatives are included, 
a sixth-order differential equation is obtained for the wavefunction and the solution by Mellin 
transform is regular in the a — > limit. It follows that inclusion of the scalar field in the 
quadratic gravity action is necessary for consistency of the quantum cosmology of the theory 
at very early times. 



The quadratic gravity theory derived from the four-dimensional heterotic string effective 
action obtained by orbifold compactification of the extra dimensions contains the coupling of 
the dilaton field to the Gauss-Bonnet term Rq B and a pseudo-scalar field to an RR term [1]. 
While the quantum cosmology of this theory can be defined on a superspace which includes 
metrics with anisotropy, the properties of renormalizability in the generalized sense and uni- 
tarity occur only when the RR term is set equal to zero. This can be achieved either by setting 
the pseudo-scalar field equal to zero without necessarily imposing any restrictions on the metric 
superspace or by confining the quantum cosmology to a minisuperspace of isotropic metrics. 
If the minisuperspace is chosen to be the space of Friedmann-Robertson- Walker metrics with 
scale factors {a(t)} and scalar fields {$(£)}, integration over the spatial coordinates gives a 
volume factor and leaves a one-dimensional action with an ad 2 term. While a boundary term 
can be added to produce an action depending only on a, a and <E>, different quantization pro- 
cedures may be used according to the choice of conjugate momentum variables. For an action 
with the second derivative of the scale factor, the conjugate momentum P„ can be introduced 
leading to a closed- form differential equation for the quantum cosmological wavefunction [2]. 
Use of the conventional momenta P a , P$ gives rise to a Wheeler-De Witt equation which 

can be expanded in powers of S - T - with the higher-order terms representing corrections to the 
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second-order Wheeler-DeWitt equation [3]. At very early times, however, the higher-order 

curvature terms will have a magnitude comparable to the Ricci scalar, so that a truncation of 

the series expansion in S - T - would no longer accurately define the Hamiltonian. Consequently, 

the sixth-order closed-form differential equation is preferable for the initial Planck era, whereas 

— * 

the series expansion of the Hamiltonian in ^-3- would be valid in the inflationary epoch. While 
this procedure might be considered to be at variance with the standard canonical quantization 
of a gravitational field theory, different methods of quantization also have been found to be 
necessary to obtain the appropriate form for the Wheeler-DeWitt equation for higher-order 
gravity [4]. It also would be consistent with the freezing of the additional degree of freedom 
during the transition from the Planck era to the inflationary epoch. For an inflationary model 
based on a slow-roll scalar potential, derivatives of the wavefunction with resepect to $ can 
be neglected [5] . This approximation can be used to obtain a limiting form of the closed- form 
differential equation near the boundary between the Planck era and the inflationary epoch, 
which is a third-order differential equation in the scale factor a. After imposing matching 
conditions on the wavefunction and its first two derivatives at the boundary, the numerical 
solution extrapolated to a = is found to be singular. An exact solution of the sixth-order 
differential equation by series expansion can be obtained by deducing the coefficients from 
a two-variable recursion relation. Techniques have recently been developed for solving two- 
variable recursion relations [6] , but the approximate solution by Mellin transform is sufficient 
to establish regularity of the wavefunction in the a — > limit. 
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1. Solution to Differential Equation for the Wavefunction near the Inflationary 
Epoch 



Given an isotropic cosmology and specializing to Friedmann-Robertson- Walker metrics, the 
low-energy effective action for the heterotic string can be reduced to a one-dimensional integral 



I = 



dt 



(6a A d + 6ad 2 + 6aK) + ^a 3 $ 2 + 6 : -^-a(d 2 + K) 

2 9a 



(1) 



with conjugate momenta [2] [3] 



Pn = 



dL d fdL\ 



da dt \dd J 



- = 6-^<S>(d 2 + K) 
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Pn = 6 a 



Pj, = 



9l 



■{d 2 + K) 



91 3 i 

— - = a 3 $ 
<9$ 



(2) 



so that the Hamiltonian is 



H = PaCl + Pa'd + P^-L 

= -6a(d 2 + K) + 6— $(d 2 + K)a + -a 3 $ 2 
9a 2 



-o^Wp 4 + ^ 



9l 
6 



P^e*PyP a - KPl 



(3) 



and the Wheeler-DeWitt equation P\& = can be transformed into 



„ 2 d 3 ^ ^ , ^ _ 2$ / d 3 * ^ 



.^e"* I n'—r-rr + 6«"— r— ) - ZvV 



4 4 

= a g 4 



4a 3 — + a 4 — T 
oa aa 2 



a#fe 



V<9a 2 <9$ 2 ~ da^d^J 

\dad<S> 3 + dad$ 2 
d 



+ 2^ 6 r5a5¥"5$3j + 4a^ 6 " 0* " 



(4) 



The correspondence between the solutions of the sixth-order equation and the pseudo- 
differential equation = can be defined as follows. Given the two operators 



A = glP^e*a 4 P a 

b = \§p^pyp a 



2a 3 



KPi 



(5) 
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the solutions to = belong to ker (A - B). If C = e - * P^e'* P$(A 2 - B 2 ) then ker C = 
ker(A 2 -B 2 )u(A 2 -B 2 y l ker(e~® P^e~^ P$). Since (A 2 -B 2 )<b = (A+B){A-B)V + {AB- 
BA)V, 

[ker(A-B) U (A-B)- 1 ker(A + B)] n ker(.AB - #4) C fcer(^ 2 - B 2 ) 
fcer(„4 - B) = ker[(A + B)(A- B)\ f)(A- By 1 coker(A + B) (6) 
D [fcer(„4 2 - B 2 ) n ker(A8 - BA)\ C\(A- By 1 coker(A + B) 

If * G fcer(.4 - then (AB - BA)^ = (A 2 - B 2 )^, so that the equivalence ker (A - 
B) = ker (A 2 - B 2 ) n (A - By 1 coker(A + B) is valid in this space. As ker (A 2 - B 2 ) = 
A^ 1 ker (A — B) n B^ 1 ker (A — B), the restriction to the subdomain ker (A — B) implies that 
ker(A 2 - B 2 )\ ker ( A - B ) = A^ 1 ker(A - B) since A^ 1 ^ = B^ 1 ^. Because 

ker(A 2 - B 2 ) = ker CD (A 2 - B 2 y l coker(e-* P^ P<s>) (7) 



the subdomain of ker C which satisfies A^ = B^ is given by 

kerCr\{A 2 -B 2 y 1 coker{e-' i 'P^e- q 'P^) r\(A-By l coker(A+B) C ker(A-B) (8) 

- ker(A-B) 

However, the spaces of functions cofcer(e~*Pj>e~ <I> Pi>) and coker(A + B) are sufficiently large 
that the intersection does not exclude a significant portion of ker C\ ker (A-B)- 



With the addition of the potential term V($), the Hamiltonian becomes 



H = -glP^e^a 4 P a + ^P 2 



9 -^p^p 2 yp a -KPi 



+ a 3 F($) (9) 



Transposing the square root expression leads to the equation 



glp-'e^a+Pa - ^P 2 - a 3 V($)) * = 
and upon squaring this equation, the new terms are 
g 2 P^e*a 4 P a - ^P 2 ] (-a 3 V($))V 



i±P^e*P 2 yp a -KP 2 a 



(-a*V(<S>))[glP^e*a 4 P a 



2a 3 



(10) 



P| )* + a 6 1/($) 2 * 



Multiplying by a single power of e Pj> gives 



g 2 a*P a a*P^e* - |p 2 e* 'aP a + ^e^Pl - | a 4 P a lp|J* 
= ^P 2 a 3 P a * - Ke^P^P 2 a ^ 



(11) 



(12) 



Since P$ is represented by the differential operator and its inverse P^ 1 by the integral 

i J d&. For a slow-roll potential 1 and given that the $ derivative of the 

potential and wave function is negligible in the semi-classical regime, 

giving rise to a third-order differential equation 

^0 + a(l + e*a 6 )^ + [(1 + c*a 6 ) - 2a«V(S)]^ " 3a 5 F($)* = (14) 

The solution to this equation can be matched with the wave function satisfying the second 
Wheeler-DeWitt equation in the inflationary epoch after specification of the boundary condi- 
tions. This equation differs from a standard equation with form [7] 

x 2 y"' + x{ax n + b + c+l)y" + [ax 2n + (ac + 0) + 7 + be] y' + (c- l)(ax 2n + I3x n + 7) - = (15) 

X 

by a single term, after setting 7 = 0, a = 0, n = 6, a = 6e*, b + c + 1 = 6, 6c = 6, 
ac + /3 = 12e* + (3 = Qe® — 12V($). It is known that the substitution w = xy' + (c — l)y leads 
to a second-order equation of the form x 2 w" + x(ax n + &)«/ + (ax 2n + /3x n + j)y = 0. 

The a — > limit of the differential equation (10) is 

a 2 d 3 ^ d 2 * d* „ 

«-J-T + "J- = (16) 



or equivalent ly 



with ^ set equal to X . If 



6 da 3 da 2 da 



d 2 d 

#0 = a 2 —: + 6a— + 6 (18) 
da z da 



the full equation is 

(H + H 1 )x = 

H lX = 6a 7 e*^ + 6a 6 (e* - 2V($)W - 18a 5 V($) / da X 
da J 

Since the solution to H xo = is xo = 4- + -3, 



(19) 



J \ a 2 a 3 / \ a 2a 2 

= 6Aa 4 (F($) - e*) - 3B(V($) + 4c*) - 18Ca 5 

(20) 



If H xi = —Hixo, then H(xo + Xi) = Hixi- With the additional terms, a series solution 
12^=0 Xn is obtained for \ with H Xn+i = —HiXi- Defining the Wronskian to be 

d d A d B B d A AB . , 

W = X01~rXV2 ~ X02~ r X01 = 3 3^ o = 6~ ^ 21 ) 

da da or da or or da or a D 

and noting that the differential operator H is «~ 4 ^a 6 ^- + 6, the solution to the equation 
HoXi = -Hixo is 

u i 7 . y, -2 (-#1x0) y. -2 (--ffixo) 

Xi = hxoi + ^2X02 + X02 / aaa xoi 777 X01 / aaa X02- 



W ^ J ^ W 

which equals 

A f B _ 2 [6Aa 4 (V{<S>)-e*)-Wa 3 {V{<S>) + 4e*) + 18Ca 5 V($)} 



Xi 



= 7j / aa^-a 
ar J cr 



AB 

a 6 



B [ A _ 2 [6Aa 4 (F($)-e*)-3£a 3 (y($)+4e*) + 18Ca 5 ^($)] 



I da -a 



a 3 J ^a? a *b 
= ^a 3 (V($) + 4e*) - -a\V{<5>) - e*) + |a 5 y($) 

after setting fci = &2 = 0, since an additional constant in the expression for each integral only 
leads to an overall shift in the value of A and B and therefore can be absorbed into those 
parameters. Then 

A B R A 

Xo + Xi = £ + ^ + ^a 3 (F($) + 4e*) - ya 4 (F($) - e*) + ^a 5 F($) (24) 

and 

ABB A "\C 

* + * x = ____ + D + i 7 j a 4 (F($) + 4e*) - ^a 5 ^*) - e*) - — a 6 F($) (25) 

Substitution of this function into the third-order differential equation gives D = 0, except at 
a = 0. Then 

4 R R 4 ^r* 

*o + *i = -^"^ + ^W*) + 4e*)a 4 - ±(V(*) - e*)a 5 - ^V(<S>)a 6 a + (26) 

At a = 0, it is possible to set D = 1 to satisfy the initial condition <]/(a = 0) = 1, although 
D = for all a would be necessary for continuity of the wavefunction. 



(23) 
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At the next order 



Hi 



n 4 n 3 QC 

~(V(^) - e*)A + Uv{3>) + 4e*)B + ^a 5 V($) 



= 6a 7 e* 



3a 



45C 



+ 6a 6 (e* - 2F($)) 

.5 



^>(d>) - e *)A + + 4e*)£ + ^a 5 F($) 

3C 



18a 5 V($) 



?L(V(*)-e*)A+^(y(*) + 4e*)B 



56 



■a 6 V($) 



(27) 



150e* + 78V($) 



= A- 



a 10 (V($)-e*) + B 



48e* - 33W$ N 



35 v v / / ■ 20 

and the solution to the equation #0X2 = —Hixi is 
X2 = -^/ 



■a 9 (^($) +4e*) 



da a 



20 



—5- [ da a 2 
a 6 J 



+ c (r $ -^))« n ^) 



A- 



-150e* + 78F($) 
35 



Cl ye 



a 10 (V($)-e*) + B 
V($) J a n 1/($) 



48e* - 33V($) 



20 



a 9 (F($) + 4e*) 



81 $ 135 \ n 



28 



= A 



(75e* - 39 1 [/($))(F($) - e*)a 10 (16e* - 11V($))(V($) + 4e*)a 9 

— /~s 



2730 



B- 



880 



(28) 



It follows that 

4 R R 4 Q(° 

X = ^ + ^ + + 4e*)a 3 - y (V(<&) - e*)a 4 + ^($)« 5 

" ^(16e $ - IWmyym + 4e*)a 9 + - 39F($))(F($) - e*)a 10 + (29) 



and 



35 



a 2V + 40^ +4 ^ 

^(lGe* - lW{*))<ym + 4 e <3 > 10 + ^(75e* - 39F($))(F($) - e*)~ n 



30030 



728 V 7 28 ' 



(30) 



This wavefunction may be extrapolated to a = a& representing the boundary between the initial 
era and the inflationary epoch. There it can be matched with the no-boundary wavefunction 



ONB 



(a) = 



by equating the derivatives with respect to the scale factor up to second order. 



(31) 



A series solution of the form c n (a — ai,) n , where a& is the value of the scale factor 
at the boundary separating the initial Planck era from the inflationary epoch, can be used to 
match the solution to the third-order differential equation with the inflationary wave function. 
A recursion relation is obtained for the coefficients of this series 



(n + 5)(n + 6)(n + 7)c n+7 + 2a h (n + 6)(n + 7)(n + 8)c n+8 

+ a 2 b (n + 7)(n + 8)(n + 9) t-n+9 



+ 



+ 



(1 - 6a^e*)(n + 6)(n + 7)c n+7 + (a 6 - (1 + age*))(n + 7)(n + 8)c n+8 

- 15a^e*(n + 5)(n + 6)c n+6 - 20ajj(n + 4)(n + 5)c n+5 e* 
15a^e' I> (n + 3)(n + 4)c n+ 4 - 6ab(n + 2)(n + 3)c n+3 e* — (n + l)(n + 2)c n+2 e* 

(1 + age* - 2a6V($))(n + 7)c n+7 + 6a£(e* - 2F($))(n + 6)c n+6 

+ 15a^(e* - 2V($))(n + 5)c n+5 + 20a3(e* - 2V($))(n + 4)c n+4 
+ 15ag(e $ - 2V($))(n + 3)c n+3 + 6a 6 (e* - 2F($))(n + 2)c n+2 

+ (e*-2F($))(n + l)c n+1 



(32) 



c n+ i + 5a 6 c n+2 + 10a 6 c n+3 + 10a 6 c n+4 + 5a b c n+5 + a 6 c n+6 



= 



with Co, c\ and c 2 given by ^m/. (o>b), ^i n f.( a b) an d f^in/. ( a *>)- These values may be deter- 
mined, for example, from the no-boundary wavefunction, whereas a&, e* and V(<&) can be 
obtained from the time-dependence of the scale factor and the heterotic string potential. 

The large-n limit of the recursion relation is a^c n+ g + 2abC n+ s + c n+ -? = with solu- 
tion c n ~ ( — at) n n and substitution of this functional dependence into the series gives 
J^r, (—ab)~ n n (a — ab) n = - — =4 r = — , which is the same divergent growth obtained ear- 

* — " \ / \ 1 l+o. (a— ah) a 



lier for the wave function as a — > 0. Given standard values for e*, V($) and a&, the recursion 
relation can be solved numerically for the higher-order coefficients, and it may be confirmed 
that divergent behaviour is again obtained for the wavefunction near a = 0. 
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Divergence at a = implies that a theory without dependence on the scalar field could not 
be physically consistent unless a lower bound is introduced for the scale factor. The existence 
of a lower bound would imply that the quantum theory must be defined on a minisuperspace 
{a(t), $(£) | a(t) > ao} so that it describes perturbations about classical solutions with spacelike 
sections of non-zero minimal radius. 

2. Solution to the Differential Equation for the Wavefunction in the Planck Era 

The change of variables w = leads to the following form for the sixth-order differential 
equation: 



6- 



d 3 * 



— a 2 w 2 'i 

6 \ da 3 dw ' " da?dw J 



Kw 3 w 



d 4 V 



d 3 $> \ 



4 4 

= a a 4 



4a 3 — +a 4 — T 
oa oa z 



<9 2 * 3* 



da?dw 2 da 2 dw ) 
a 9i w I Q a Q w 2 W dadw 3 ) 



2 2 

2^4 ^ dadw ' dw 



+ — + 3u> 



d 2 ^ 



,d 3 ^ 



w 



(33) 



i 3 fm_ 16 Q 2 ^ 

4a 6 V <9io <9u> 2 



d-u) 2 9i« 3 / 

<9io 3 Sir 4 



d 5 * 



+ 115io 2 — + 75ty 3 -— r + 16w 4 — + it) 



<9io 5 <9u> 6 / 



Given that the Mellin transform has the property A4[f^; s] = (— l) n f*(s — n), \t*(s,a) 

satisfies the equation 



-a 29 j(s + l) 



' d 3 **(s + l,a) <PV*(8 + l,a) 



da 



da 2 



4 4 

a 4 



4a' 



,d**(s,a) 4 d 2 **(s,a 



da 



da 2 



+ if(s + l)(s + 2) 
a#l(s + l) 2 (s + 2) 



d 2 **(s + 2,a) 
da 2 

dfr*(s + l,a) 
da 



3 2 



a( g + l) ^±M + (g + l)3^ + l,a) 



(34) 



+ — r[s 6 + lis 5 + 50s 4 + 125s 3 + 205s 2 + 242s s + 152]** (s + 2, a) 
4a° 

The difference-differential equation is 
L 1 (s,a)E^*(s,a) + L 2 (s,a)E 1 ^*(s,a) + L 3 (s,a)^f*(s,a) = 

LAs,a) = -K(s + l)(s + 2)Dl + -*-r[s 6 + lis 5 + 50s 4 + 125s 3 + 205s 2 + 242s + 152] 

4a b 

L 2 (s, a) = a 2 |(s + l)(aD 3 + 6D 2 2 ) + ^ 2 (a(s + 1)D 2 + (s + l) 3 ) 

+ agl(s + l) 2 (s + 2)D 2 



L 3 (s,a) = a?gt[Aa 6 D 2 + a 4 ^ 1 
£7i**(s,a) = **(s + l,o) 



D 2 **(s,a) = ^*(s,a) 



which has the asymptotic form 
1 



-^s 6 u s+2 + agj(s + l) 2 (s + 2)D 2 u s+1 + -gl(s + l) 2 u s+1 ~ 



(36) 



with solution 



Us+2 
U s +1 



(4a 6 ) 



u s+2 ~ (-!)' 



(4a 6 )' , , 
YW a9 * D2+ 2 9 * 



U 2 



The equation containing u 2 is 



(37) 



-KD%u 2 + ^ 152 ^2 + a 9 ^(aDl + QD 2 2 )u l + ^gl ( -■//>, ! I ) a, 



3 o /5 

6 ' —-"-^ ' 2^* V 3 

+ a 4 ^(4a 3 D 2 + a 4 Dl)M = 



(38) 



so that 



u 2 



KDn 



38' 



-l 



+ 6D%)ui 



3 2 

of4 



raJ D 2 + 1 ) «i + a 4 gj(4a 3 D 2 + a 4 Df)u j 

(39) 



For large s, the Mellin transform of the wavefunction is 



(4a 



6\s-2 



ag\D 2 



3 2 
2^4 



s-2 



-1 



.51 



a^(aDi + 6Di)**(l,o) 



+ \ (^aD 2 + 1 ) * (LrO + r/ L r/}(4r/ :i i9 2 + « '^)**((J. «; 



Application of the inverse transform yields the wavefunction 



(40) 



ty(w,a) ~ 



(4a 



6\s-2 



ag\D 2 + -g\ 



s-2 



w s ds 



KOI - 5 



a 9 l{aDl + 6Di)**(l, a) + ^ QaL> 2 I 1 ) 'I' I • " 



+ a 4 <7 4 (4a J D 2 + a 4 ^)**(0,a) 



(41) 



After shifting the variable s to s — 2 and choosing the new contour to be the imaginary axis, 
the integral has the form 



w 



—2 r°° 



2tt 



T(-iy) 3 (-4a 6 ) iy 



y sink ny 



IT 



ag 4 D 2 + -0 4 



iy 



W 



-iy 



dy 



(42) 



or equivalent ly 



w -2 roo 



16vr 4 



/oo 
T(-iy) 3 y 3 e R+4M-4a 6 )+M«9l £>2+f sl)-^ *>)]]v dy 
-OO 

-2 poo 

I T(-iy) 3 y 3 e ^+ i [ ln (- 4a6 )+ ln ( a 9 2 4 D 2 +^ g 2 4 )-in w)]]y dy 

J — oo 

3W f°° 3 3 [- 7T+ i[l n (-4a 6 )+ln(aglD 2 + ^g 2 4 )-ln w)]]y . 

16vr 4 7^ 1 y; ^ ^ 

/oo 
r(-zy) 3 y 3 e [-3vr+i[Zn(-4a 6 )+«n(a 9 |D 2 + |< ? l)-;n tu)]]j, rfy 
-oo 



3u> 

167T 4 

-2 />oo 



16vr 4 /_ r 



The function 



- se< s 



Y{s) 11 V n 

v ' n=l 



oo 



has no poles or essential singularities in the finite plane, and it decays at a factorial rate 
if Re s > 0. Consider the contour obtained by rotating the upper and lower halves of the 
imaginary axis towards the real axis. Since the integrand vanishes on the circular contours C\ 
and C2 , the integral over the imaginary axis equals the integral over a contour which traverses 
the positive real axis in opposite directions and therefore equals zero. 

Given that equation (40) defines a close approximation to \t*(s,a) for \s\ > N , a more 
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accurate formula for the wavefunction is 



1 f io ° (4a 6 ) 3 - 2 



2) ; 



s-2 



w s ds 



2 a 6 



a 9 j-(aD 3 2 + 61>f )**(!, a) + \g\ (^aD 2 + 1 ) * ( \ .a) 



a A gi(4a i D 2 + a A D'^*(0,a) 



iNn 



/ ^*(s,a)w s ds 



1 

2vri 



iNn 



{ 1 n* - 2) 3 



Cig 4 D 2 +2^4 



s-2 



w s ds ■ 



2 a 6 



a 9 l{aDl + 6D 2 2 )**(1, a) + \g\ QaD 2 + 1 ) **(!, a) 



a 4 5 4 (4a 3 L> 2 + a 4 !^) 



4n2 



i**(0,a)} 



1 

1 

2vri 



iJV 

^*(s,a)w~ s ds 

■iN 

iN ° , 11s (4a 6 r 2 
- Wo r (* " 2) 3 



3 

a5 4 A> +2^4 



s-2 



io s ds ■ 



KD 



38 



2 Tfi 



-i 



+ 6I$tf*(l, a) + ^-aD 2 + 1 ) tt* (1, a) 
+ a 4 </ 4 (4a 3 J D 2 + a 4 L> 2 )**(0, a) 



(45) 



because the first integral vanishes. The values of 5 f *(s, a) on the imaginary axis have not been 
given, but they can be defined by analytic continuation from the values on the real axis, which 
may be deduced from {\I/*(s,a)| — 2 < s < 2} through the recursion relation. The boundary 
conditions for the wavefunction can be transformed to conditions on the Mellin transform. 



{s,a b ) 
(s,a b ) 



/» OG f'OQ 

Jo Jo 



inf. 



da 
da 2 



da 



(s,a b ) 



(46) 



inf. 

da 2 



(s,a b ) 



Given ^* n f (s,a b ) for arbitrary values of s, the the higher derivatives ^*(s,a b ), — 2 < s < 2, 
can be obtained from the difference-differential equation (34) and ^*(s,a), —2<s<2, can 
be extrapolated to a = 0. 
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The a — > limit of the integral (41) is given by 

—2 /-oo 



w 



2vr 4 /_ 



r(-iy) 3 (-4a 6 ) ij/ y 3 sma 3 7ry 



ag 4 D 2 + -5 4 



1 a 6 



J±{aD\ + 6P> 2 )**(l,a) + ^ QaP> 2 + 1 ) tt*(l,a) 



a 4 # 4 (4a 3 P> 2 + a 4 £>£)**(0,a 



4 7-i2\,t,* , 



ID 



2 /-co 



OC 



2tt 4 



r(-iy) 3 (-4)^<5(y)y 3 sm/i 3 7ry 



3 2 

o#4 



ly 



w iy dy 



Urn,, ^D 2 2 \ a 9 4{aDl + 6P>f)**(l, a) + ^ ( -r/i9 2 f I j ¥' (La) 



3 
2' 



+ a 4 < 7 4 (4a 3 A 2 + a 4 P 2 )**(0,a) 



(47) 



and again the integral over y vanishes. Then 



lim a ^o^(w,a) = Um a ^o / ty*(s,a)w s ds 

2™ J-iN 



— lim a ^Q / (-1 



(4a 



6\s-2 



r(s-2) 3 



3 
T 



kd i - 3 4 




a . 





a0 4 D 2 + 7^4 

ci . o o . 3/5 



s-2 



<7 4 2 ( -{aDl + 6D 2 ) + ~- I |aL>2 + 1 ) ) **(l,a) 



+ a 4 a 4 (4a 3 Z> 2 + 



a 4 D 2 )vT(0,a)} 



lirria^o / a)w~ s ds 

2ni J-iN 



(48) 

because the second integral vanishes by the regularity of the integrand in the a — > limit as 
any divergences in \I>*(0, a) and a) could be absorbed by the factor of a 6 multiplying the 

operator D^ 2 . This implies that the quantum theory would be consistent in the a — > limit 
when the scalar field degree of freedom is included. 

Since the new terms in the sixth-order differential equation are given by 

1 



e -*P* e -*P* ( gjP^e*a*P a - ^ P 2 . ) (-a 3 y(d>))* 



2a 3 ' 



(49) 



P|)tf + e^P^^a 6 ^) 2 * 



regularity of the integrand in the expression for ty(w,a) and therefore consistency of the 
quantum cosmology of the theory in the a — > limit will be maintained in the presence 
of a potential V(<&). 
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